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ABSTRACT 

A duality  theorem  of  Heathcote  exhibiting  a relationship 
between  first  passage  times  of  the  queue  length  In  the  GI/M/1 
queue  and  the  busy  period  of  Its  dual  M/G/l  queue  Is  generalized 
to  the  phase  type  queues  GI/PH/1  and  PH/G/1.  The  phase  type 
distributions  Include  a number  of  well-known  models  such  as 
generalized  Erlang  and  hyper-exponential  as  special  cases  and 
form  a versatile  class  with  a number  of  Interesting  closure 
properties. 

KEY  WORDS 

Queueing  theory,  duality,  phase  type  distributions 


1.  Introduction 


The  problem,  discussed  In  this  paper,  arose  in  the  context 
of  a larger  investigation  into  the  nature  of  the  busy  period  of 
subclasses  of  the  GI/G/1  queues,  which  are  denoted  by  the 
symbols  PH/G/1  and  GI/PH/1  and  which  are  defined  below.  We 
attempted  i.a.  to  generalize  the  well-known  duality  result  of 
Heathcote  [2],  which  relates  the  limit  distribution  of  the  times 
between  points  of  increase  for  the  maximum  queue  length  process 
in  the  GI/M/1  queue  to  the  distribution  of  the  busy  period  of 
the  dual  M/G/l  queue.  This  result  implies  in  particular  that 
in  an  unstable  GI/M/1  queue  the  maximum  queue  length  process 
grows  approximately  like  the  counting  process  of  a renewal 
process . 

The  duality  theorem  of  Heathcote  carries  over  to  the  GI/PH/1 
and  PH/G/1  queues,  but  the  technical  difficulties  of  the  proofs 
are  substantially  greater  than  in  [2].  Several  deeper  properties 
of  the  matrix  transforms  used  in  the  study  of  these  queues  are 
needed,  and  the  methods  of  the  present  paper  may  be  of 
independent  interest.  Our  main  theorem  implies  in  particular 
that  the  maximum  queue  length  process  in  an  unstable  GI/PH/1 
queue  grows  a pprox ima tely  like  the  counting  process  of  an 
appropriately  defined  Markov  renewal  process. 

2.  Phase  Type  Distributions  and  Phase  Type  Renewal  Processes 
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distribution  on  [0,«)  which  is  obtainable  as  the  distribution 
of  the  time  till  absorption  in  a continuous-time  finite  state 
space  Markov  chain  with  a single  absorbing  state  into  which 
absorption  is  certain.  The  class  of  such  distributions  includes 
a number  of  well-known  particular  cases  such  as  generalized 
Erlang  and  hyper-exponential  ( 1 . e . , a mixture  of  a finite  number 
of  exponentials)  distributions  and  due  to  its  Interesting 
closure  properties  [4]  constitutes  a versatile  class  with 
properties  especially  useful  in  the  algorithmic  solution  of 
several  queueing  models  (c.f.  references  in  [5]  and  [8]). 

To  be  specific,  consider  a Markov  chain  with  state  space 
{ 1 , . . . ,m ,m+ 1 } , initial  probability  vector  (a,am+^)  and 
infinitesimal  generator 


where  o«(aj,...,a|n),  T is  a non-singular  mxm  matrix  with  T^<0 
and  T^zO  for  i ^ j , and  T°>0  is  an  m-vector  satisfying  Te+T°=0 
with  e ' = ( 1 , . . . , 1 ) . For  such  a Markov  chain  the  probability 
distribution  of  X,  the  time  till  absorption  in  m+1  , is  given  by 
the  c . d . f . , 

H ( x ) = 1 -oexp(Tx)e,  x>0  (2.1) 

Def 1 ni ti on  2.2:  Any  probability  distribution  H(*)  obtained  as 
above  Is  called  a Phase  Type  Distribution  ( nM - d 1 str i but  1 on ) , and 
the  pair  (a,T)  is  a representation  of  H(*). 

Remark ; To  avoid  uninteresting  complications.  In  the  sequel  we 
shall  assume  c»m+,=0  so  that  H(*)  does  not  have  an  atom  at  0. 
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Example  2.3:  The  generalized  Erlang  distribution  which  Is  the 
convolution  of  m independent  exponential  distributions  with 
parameters,  say,  u.,...,p  respectively  has  representation 

a* ( 1 ,0  ....  ,0 ) 


-U 


1 


1 


0 


-u. 


0 
0 

" l'm 

Example  2.4:  The  hy per- exponent  1 a 1 distribution  which  is  the 
mixture  of  m exponentials  with  parameters,  say,  vij,...,vim  has 
a representation  T*d1 ag ( -p j , . . . ,-p  ) with  the  components  of  a 
giving  the  respective  mixture  ratios. 

b.  Renewal  Processes  of  Phase  Type: 

Each  time  the  Markov  chain  Q becomes  absorbed  in  the  state 
(m+1),  restart  it  by  performing  a multinomial  trial  with  possible 
outcomes  l,...,m  and  probabilities  c^,...,am  to  pick  a new 
"initial  state".  Considering  each  absorption  into  the  state  m+1 
as  a renewal,  we  obtain  a renewal  process  for  which  the  time 
between  any  two  successive  renewals  has  c.d.f.  H(«)»  the  PH- 
distrlbution  given  by  (2.1).  Such  a renewal  process  is  called  a 
Renewal  Process  of  Phase  Type  (PH-renewal  process)  (Neuts  [5]). 

The  above  procedure  also  constructively  defines  a new  Markov 
chain  with  state-space  {l,...,m),  initial  probability  vector  a 
and  Infinitesimal  generator  Q*=T+T°A°,  where  A°adiag(aj  ,. . . ,a  ) 
and  T°*(T° , . . . ,T° ) . This  Markov  chain  describes  the  "phase"  of 
the  system  and  is  of  considerable  importance.  In  [4]  it  is 
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shown  that  one  may,  without  loss  of  generality,  assume  that  the 
representation  (a.T)  of  H(*)  Is  so  chosen  as  to  make  Q* 
Irreducible,  and  we  shall  henceforth  assume  that  this  is  Indeed 
the  case. 

We  let  e denote  the  stationary  probability  vector  of  the 

Markov  chain  Q* , l.e.,  the  unique  (strictly  positive)  vector 

satisfying  oQ*=0,  ee=l.  It  may  be  easily  verified  that  o^-XaT'1 

where  X~^=-aT"^e  Is  the  mean  of  H ( • ) . 

The  m*m  matrices  P(v,t),  v>0,  t>0,  defined  in  [5]  are  such 

that  the  entry  pjr<”  ,t)  is  the  conditional  probability,  given 

that  the  initial  phase  Is  j,  that  at  time  t+,  the  Q*-cha1n  Is  In 

state  j'  and  that  v renewals  have  occurred  In  ( 0 , t ] . It  Is 

known  [5]  that  these  have  generating  function 
def  “ 

P(z,t)  = E zvP(v,t)  = exp[(T  + zT°A°  )t]  , f z ( ^ I » t^O. 

v = 0 

We  also  recall  [6]  that,  under  the  assumption  Q*  is  Irreducible, 
the  matrices  P(v,t),  v>l,  t>0,  are  all  strictly  positive. 

Further,  it  may  be  easily  seen  that  P^(0,t)>0  for  all  t>0, 

1 < 1 <m . 

s s 

3.  The  Phase  Type  Queues 
a.  The  GI/PH/1  Queue: 

Consider  a GI/G/1  queue  In  which  the  service  time  c.d.f. 

H ( • ) is  of  phase  type  with  representation  (a,T)  and  where  the 
Inter-arrival  times  are  1 . 1 . d . , with  a non-degenerate  probability 
distribution  F(*).  Such  a model,  denoted  by  GI/PH/1,  has  been 
discussed  In  detail  by  Neuts  [6],  and  for  later  use  we  quote  the 
following  results. 


Let  t denote  the  epoch  of  the  n-th  arrival,  c the  size  of 
n r n 

the  system  at  tn+  and  the  phase  of  the  service  at  rn  + . Then 
^*n*^n,Tn~Tn-l  ^ 5 where  Tq"t_j"0,  defines  a semi-Markov 

sequence  with  state  space  { 1 ,2  ....  }x{  1 ....  >}*[(),•)  and  transition 
matrix  Q ( x ) given  by 


Q ( x ) 


Bq( x ) Aq(x  ) ...  ... 

B j ( x ) A ^ ( x ) A q ( x ) 

B2 ( x ) A2(x)  A ^ ( x ) AQ(x) 


. x>0. 


where , 

x 

A (x)-/P(n,t)dF(t),  n>0 , x>0, 

‘ A 


" X 

B ( x ) ■ l /P ( v , t ) dF ( t ) A°  ° , n>0 , x>0, 

n v-n+1  0 


where  A°°  Is  the  mxm  matrix  each  of  whose  rows  Is  a. 

Since  F(*)  Is  non-degenerate,  we  have  (from  the  strict 

positivity  of  the  matrices  P(v,t)  for  v>l,  t>0  and  the  positivity 

of  the  entries  P^(0,t)  for  t>0  mentioned  at  the  end  of  Section  2) 

de  f - 

that  the  matrices  An  ■ An(«),  n>l,  are  all  strictly  positive; 
further,  the  matrix  A0d-fAQ(«>)  Is  such  that  each  of  Its  diagonal 
entries  Is  positive.  This  entails  that  the  embedded  Markov  chain 
with  transition  probability  matrix  Q(«°)  - with  possibly  some  of 
the  states  (l,j)  removed  - Is  Irreducible  and  aperiodic.  We 
shall  conclude  our  Introduction  to  the  GI/PH/1  queue  by  recalling 
[6]  that  this  queue  Is  stable  Iff  p-63>l,  where  e.  Is  (also)  the 

def  " 

Invariant  probability  vector  of  A ■ T.  An,  and  I nAne. 
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b . The  PH/G/1  Queue  : 

By  PH/G/1  we  denote  a Gl/G/1  queue  in  which  the  arrival 
process  Is  a PH-renewal  process.  Note  that  this  queue  Is  the 
dual  of  the  GI/PH/1  queue  defined  above  In  that  It  may  be 
obtained  by  reversing  the  roles  of  1 nter-a rr 1 va 1 times  and 
service  times  In  the  latter.  The  PH/G/1  queue  Is  a sub-class 
of  the  more  general  N/G/l  queues  studied  by  V.  Ramaswami  [8], 
and  we  shall  state  below  some  basic  properties  concerning  the 
PH/G/1  queue  by  particularizing  results  obtained  in  [8]. 

Let  us  assume  that  the  arrival  process  is  the  PH-renewal 
process  with  representation  (a,T)  and  that  the  service  time 

c. d.f.  is  given  by  F ( • ) . Defining  tn  to  be  the  epoch  of  the  n-th 
departure  (r^O),  and  Cn  and  Jn  to  be  respectively  the  queue 
length  (i.e.,  the  number  of  customers  In  the  system)  and  the 
phase  of  the  arrival  process  at  1^+,  It  is  easily  seen  that 


. * A 

n*?n  + r'n): 

n >0 } 

is  a semi-Markov 

sequence  with  state 

space 

{0,1  .... > x { 1 . . ,m ) x [ 0 ,«) 

and  trans 1 tion  matri x 

X 

o 

L’° 

Cj  (x) 

C 2 ( x ) .. 

. 

A0(x) 

A 1 ( x ) 

X 

C\J 

• • • • 

0 ( X ) - 

0 

Aq(x) 

A] ( x ) .. 

• • • • 

, x>0 , 

0 

0 

%( * ) 

• • • 

• • • 

*— J 

where 

An( • ) , n>0 , 

are  as 

def 1 ned 

earlier 

in  (a ) , 

and 

C(x)«/exp(Ty)‘T°A°*A(x-y)dy,  n>0. 

" 0 

It  Is  easily  seen  that  the  Markov  chain  defined  by  Q(«)  Is 
irreducible  and  aperiodic.  We  recall  from  [8]  that  the  PH/G/1 


I 

I 

queue  is  stable  iff  p"8{J<l.  Further,  p»Xu,  where  u is  the  mean 
of  F(*),  is  the  traffic  intensity  of  the  PH/G/1  queue. 

Basic  to  the  discussion  of  the  PH/G/1  queue  are  the  first 

I 

passage  times  from  the  set  of  states  W ( 1, j ) : l<j<m)  to  the 
set  of  states  0"{(0,j):  1 < j < m ) which  are  governed  by  the  matrix 

G(x),  x>0,  whose  (j.j')-th  entry  is  the  probability  that  starting 


(3.4) 
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(iv)  If  p<1,  then  G Is  stochastic.  If  p>l,  then  at  least  one 
component  of  Ge  is  less  than  one. 

(v)  For  any  s>0,  G(s)  is  strictly  positive  and  is  the  monotone 

oo 

limit  of  the  non-decreasing  sequence  {G  ( s ) } defined  by 

0 


Gq(s)  = 0 

CD 

Gn  + 1(s)  = E Ay(s)G^(s),  n> 0 
v = 0 


(3.5) 


Proof : All  the  results  above  follow  from  Theorem  2.2.11  In  [8] 

which  in  turn  is  proven  by  Neuts  [7]  in  a more  general  context 
in  the  analysis  of  Markov  Renewal  Processes  with  Q(*)  of  the 
form  given  above. 

Remarks : 

1.  In  view  of  the  structure  of  Q(*).  we  have  •*.  i * a t £(•)  also 
describes  first  passages  from  the  set  i + 1 = { ( i + 1 , j ) : l<j<m) 
to  the  set  i={(i,j):  l<j<m)  for  any  i>0. 

2.  Note  that  the  i-th  entry  of  G(s)e  is  the  Laplace-Stieltjes 
transform  of  the  busy  period  starting  with  one  customer  and  in 
phase  i . 

3.  Equation  (3.3)  generalizes  Takacs'  equation  for  the  M/G/l 
queue . 

4.  First  Passage  Times  in  a GI/PH/1  Queue 

Consider  the  GI/PH/1  queue  defined  in  Section  3a.  By 
level  ji  we  denote  the  set  of  states  {( n , 1 ),...,( n ,m )} . The 
principal  objects  of  study  of  this  paper  are  the  first  passage 
times  from  n to  n + k in  the  GI/PH/1  queue;  these  will  be  the 
subject  matter  of  this  section.  In  the  sequel,  we  let  for  x>0, 


i 
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n,k>l,  Dn  n+k(x)  denote  the  mxm  matrix  whose  (i.j)-th  entry 
Dn  n>|t(1»j;x)  Is  the  conditional  probability  that  the  process 
enters  n + k for  the  first  time  at  or  before  x by  visiting  (n+k,j) 
given  that  it  starts  at  level  n and  in  phase  i.  We  also  define 
the  Lapl ace-Sti el tjes  transform 

Dn,nU(s)  ‘ 'e'S’‘d°n,n+k(x)-  SS°- 


| 

i 

I 


| 

! 


Proposition  4.1:  For  n>l,  k>2,  we  have 

D =D  *rj  * *D 

un,n+k  n,n  + l n+1  ,n  + 2 n+k-l.n+k* 

and 

^n,n+k  ^n,n+l  ^n  + 1 ,n  + 2 ^n+k-l.n+k’ 

where  * in  (4.2)  denotes  matrix  convolution. 


(4.2) 

(4.3) 


Proof:  The  transition  from  n to  n+k  can  occur  only  along  a path 

of  first  passages  from  n to  n + 1 , from  n + k- 1 to  n + k . Given 

any  sequence  j i , . - . » j k ^c{l,...tm)  denoting  the  phases  at  the 
epochs  of  such  first  entrances  to  n + 1 , . . . ,n  + k-l  , the  duration 
between  these  first  entrance  times  are  conditionally  independent. 
This  is  an  immediate  consequence  of  the  Markov  property  at 
transition  epochs  for  the  Markov  Renewal  Process  defined  by 
Q(*).  Hence  (4.2).  (4.3)  follows  immediately. 

Theorem  4.4:  For  n>l , s>0, 

Dn,n*l<s)  = #0<s>*nj/v<s>0n-vH.n+l<s>  + jnflv‘s>A°"Dl.n+l<s>  <<-5> 

(where  the  second  term  In  the  right  side  is  taken  to  be  0 when 
n*  1 ) . 

Proof : Let  N denote  the  number  of  departures  in  the  first 


L 
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inter-arrival  interval.  Also  let  us  denote  by  S the  event 
that  the  process  enters  n + 1 at  or  before  time  x and  that  the 
phase  at  which  it  enters  n+1  Is  j.  Then  by  a simple  probabilistic 
argument,  we  have 

P(S,N=v|c0=n,J0=1)= 


/P44(0,u)dF(u), 


if  v=0 


x m 


/ ^Pik(v,u)Dn_v+i  >n  + i U.J  ;x-u)dF(u) , if  l<v<n-l 


x ® m 


/IE  P 4 . ( v , u ) I o.D.  +1(k,j;x-u)dF(u),  If  v=n 
0 v=n  t-1  k=  1 k ,n  1 


whence 


6n.n+l(i^;x>' 


x n-1  x m 

/P.  . (0,u)dF(u)  + I / I Pik(v,u)D  , , (k,j  ;x-u)dF(u) 

0 v=l  0 k = 1 1K  n-v+i  ,n+ 1 


» x m m 

+ I / I P1t(v.u)  ^ Cl|c^l  n + 1 

v=n  0 t=l  11  k = 1 k 1 1 


( k , j ; x - u ) d F ( u ) . 


Computing  the  Lapl ace-Stl el t j es  transforms  of  these  and  putting 
the  result  in  matrix  notation,  wo  get  (4.5). 


From  the  1 rreduc i bi 1 i ty  of  Q ( - ) it  is  clear  that  the  matrix 
Dn  n+](0)  should  be  stochastic.  That  this  Is  indeed  the  case  is 
verified  below  by  directly  using  (4.5). 


Lemma 


4.6:  1S  stochastic 


Proof:  By  (4.5), 


D1 f2(0)  = A0+(A-A0)AooD1 >2(0) 


(4.7) 


AooD1>2(0)e  = AooA0e+Aoo(A-A0)AooD1 >2(0)e. 
Letting  A°°D1  2(0)e=u£,  this  yields 

u = (“Ao-)+u_ (“Ao-)-’ 


and  since  aA0e>0,  u=l.  Thus 


A° 0 D-j  >2(0)e=e 

Now,  multiplying  (4.7)  by  e and  using  (4.8),  we  get 


(4.8) 


D1>2(0)e  = AQe+(A-A0)e  = e. 

Theorem  4.9:  For  all  n,k>l,  Dn  n+(((0)  is  stochastic.  Also  for 
n>2,  k>l,  Dn  n+k(0)  is  strictly  positive. 

Proof : In  view  of  (4.3)  it  suffices  to  prove  the  results  for 
k=l.  Let  n>l  and  assume  as  inductive  hypothesis  that 

2(0),...,Dn  n+1^  are  a11  stochastic.  We  now  show  that 

Dn+1  n+2^  stochastic. 

Using  (4.5)  and  (4.3)  it  is  easily  seen  that 


EI-HJ0n+1,B+2(°)  * a„, 


(4.10) 


where 


v=Z  v = 0 

(with  the  second  term  on  the  right  being  0 when  n=l). 

Now  M>A^>>0,  and  Me=e-Age<<e.  Thus  the  strictly  positive 
matrix  M has  spectral  radius  n less  than  1.  That 


O-n)*Dn  + ltn  + 2(0)e  - 

“ ^(e-Me)  * (1-n) 

and  since  n<l,  £°n+i  Thus  Dn+1  n+2^-"-’  and  the 

proof  is  complete  by  mathematical  induction. 

Before  concluding  this  section  we  list  some  simple  but 
useful  results  as 

Proposition  4.11  : Let 

oo 

$ ( s ) **  /e~sxdF  ( x ) , s>0 . 

0 

For  all  s >0 , n , k>  1 , 

a)  Dn,nH(s)^ 

b)  Dn>n+1(s)e<t(s)e 

c)  Dnn+k(s)e<{*(s))ke 

Proof:  (a)  is  Immediate  from  the  definition  of  D_ 

n ,n  + k 

(b)  is  got  by  applying  (a)  in  (4.5).  Finally  (c)  Is  got  by 
applying  (b)  in  (4.3). 

5.  The  Duality  Theorem 

In  this  section  we  prove  the  following  duality  theorem 


i 


I 


which  generalizes  the. result  of  Heathcote  [2]  to  the  phase  type 
queues,  discussed  in  Section  3. 
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Theorem  5.1  (Duality  Theorem): 


Let  $>0  be  fixed.  As  n-*®,  the  matrices  Dn  n + ^(s)  describing 
the  first  passage  times  from  to  n+1  in  the  GI/PH/1  queue  converge 
to  the  matrix  G(s),  where  G(s)  is  the  unique  non-negative 
solution  of  the  non-linear  matrix  functional  equation 


G ( s ) = r A ( s )GV ( s ) (5.2) 

v=0 

and  which  describes  the  first  passage  times  from  1_  to  0 in  the 
dual  PH/G/1  queue.  Further  the  limit  G(0+)  is  stochastic  iff 
p<1,  i.e.,  iff  the  GI/PH/1  queue  is  unstable. 


The  lengthy  proof  of  Theorem  5.1  is  accomplished  in  several 

stages.  In  the  ensuing  discussion  s>0  is  assumed  to  be  fixed. 

Also  we  let  lim  Dn  n + -j(s)  and  1 im  D ,(s)  denote  the  matrices 

n -*■<*>  * n->®  * 

whose  respective  (i,j)-th  entries  are  1 im  Dn  ,(i,j;s)  and 

n-*«  ’ 

— Dn  n+l^’^;s^  From  Proposition  4.11(b)  it  is  clear  that 
n-*-«  * 

these  matrices  are  strictly  substochastic;  further, 


0.<HS!Dn>n+,(S)<TO  D +,(s)  (5.3) 


Lemma  5.4:  Define 


M-|  ( s ) = Aq(s) 

n n+1 

Mn+l(s)=%(s>+  1 n 

n 1 U v=l  v k=n-v+2  k 


where  the  matrix  product  is  formed  by  taking  the  terms  in 
increasing  order  of  the  index  k.  Then 
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(i)  Mn(s)  strictly  substochastic  for  every  n>l, 

( 1 i ) Mn(s  ) Is  entrywise  non-decreasing  In  n, 

(111)  Lim  Mn(s)eG(s),  where  G(s)  is  the  unique  solution  of 

n ►<*> 

(5.2);  the  limit  Is  taken  entrywise. 

Proof : It  Is  obvious  that  M ^ ( s ) is  strictly  substochastic. 

OD 

Now,  since  E A (s)e<<e,  there  exist  constants  w(s)>0  and  6(s)>0 

\>=  0 v 

ou 

such  that  w(s)+6(s)<l,  E A (s)e<w(s)e  and  An( s )e<5 ( $ )e . Now, 

v=l  v ~ u 

,-l 


M2(s)=[I-A] (s)]  A0(s)=  E A^(s)A0(s), 

v = 0 


whence 


In  other  words,  M2(s)  is  also  strictly  substochastic,  and  clearly 
M2(s)>M1(s). 

Now,  assume  as  inductive  hypotheses  that  ( s ) + ^ ( s ) and 

Ml<+1  (s  )e<<e » for  l<k<n.  We  have, 

Mn  + 2*s  ^ = ^ I_ {A1 (s  )+Az(s )Mn  + 1 (s)  + . . -+An  + 1 (s)M2(s)...Mn  + 1(s)}]'1AQ(s) 
>[I-(A1(s)+A2(s)Hn(s)+...+An(s)M2(s).. .Mn(s)}]_1A0(s) 

*Mn+ l ( s ) • 

The  inequality  above  is  obtained  by  writing  the  inverse  as  a 
power  series,  dropping  the  term  An+1 (s  )M2 ( s ) . . .Mn  + 1 ( s ) , and 
using  the  Induction  hypothesis  Mk ( s )<Mk  + ^ ( s ) , l<k<n.  Furthermore, 

OD 

Mn  + 2(s)e*  E (A1(s)+A2(s)Mn  + 1(s)  + ..  . +Ap  + 1 ( s )M? ( s ) . . .M  ^ (s ) )vA0(s)e 

vE  0 

s ««« 

by  using  the  substochastic  of  ( s ) , . . . ,M  + j ( s ) assumed  as 


Inductive  hypothesis. 


Having  proved  statements  (1)  and  (11),  at  this  point  these 


also  Imply  the  existence  of  11m  M_(s)=M(s ) . So  all  that  we  need 

n>» 


to  prove  Is  that  M(s)=G(s).  We  have. 


M„*,(s)  » A0(s)+A1(s)Mn+1(s)+...+AM(s)H9(s)...M„xl(s) 


n 2 


* EnAv(s)Mn+l(s) 
v=  0 


< Z A ( s )MV ( s ) 
v=0 


Letting  n>“>  this  yields 


M ( s ) < z A ( s )MV ( s ) 


(5.5) 


Also 


yields  upon  letting  m-*®,  that 


M(s)  > Z A (s)Mv(s) , for  all  n>0. 

v=0  v 


Since  the  above  inequality  holds  for  every  n. 


M ( s ) > Z A ( s )MV ( s ) 


(5.6) 


Now,  by  (5.5)  and  (5.6)  , 


M(s)  = Z A (s)Mv(s), 
vs0  v 


and  It  follows  that  M(s)=G(s)  by  appealing  to  the  uniqueness 


of  the  solution  to  (5.2)  stated  In  Theorem  3.2(H). 


Lemma  5.7: 


H*  V„,1<S>  5 G(s) 


(5.8) 


1 
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Proof : It  suffices  to  show  that  Dn  n+1 (s)>Mn(s)  for  all  n>l 

where  { Hn ( s ) ) are  as  In  Lemma  5.4.  Clearly,  Dj  2(s)>Mj(s). 

Assume  as  inductive  hypothesis  that  Dn  n + i ( s ) £Mn ( s ) for 
n»l , . . . ,k-l  . Then 

Dk,k+l*s)mtI'{AlM+A2*s*°k-l,kM  + ,',+Ak-l(!0D2,k(0}J'lAo(0 

> C I - ( A 1 (s)+A2(s)M|{_1  (s  ) + . . .+Ak_j  (s)M2(s)  . . .Mk  .J  (s)  > D " 1 A0  ( s ) 
=Mk(s), 

where  the  inequality  is  obtained  by  using  the  induction  hypothesis 
and  (4.3)  In  the  series  expansion  of  the  inverse.  The  proof  is 
now  complete  by  mathematical  induction. 


Lemma  5.9:  Let 


Then 


R<s>  ■ "» 


n->“ 


R(  s ) < z A (s)Rv(s) 
v=0 

Proof : For  sufficiently  large  n and  fixed  N<n-1  , we  have, 

Dn,n*l(s)’A0<s)*  s'Av(s)Vv*l,n*l(s)t  1 VS,A°°D1  ,n*I (s)  by  (4'5> 

v=l  v=n 

-Ao(s)*  " Av(OD„.vt)  ,„.,(*)♦  V A <s>Vv*l.n»l(s) 

v=l  v=N+ 1 


* f _*vCs>a*-o1.b+1(») 


v = n 


and  letting  n-*-«. 


H n-1 

R(s)<An(s)+  E A (s)Rv(s)+11m  Z A (s)D, 


for , 


I4  n v a / n \ o / • i i hi  u n l 3 

v'.  _ _ i, . i v n-v+1  ,n+l 

v=  1 n-*-»  v=N+l 


11m  Z Av(s)A°°D1  n+1(s)=0. 
n-*®  v=n  * 


(s). 


Now  letting  N-*®,  we  get 


since. 


R(s)  < Z A (s)Rv(s)  , 
v = 0 


Tim  TTm 

N-*-»  n-*-°° 


n- 1 


l A 
v = N + l v 


(s)D 


n-v  + 1 ,n  + l 


(s)  = 0 


as  is  seen  from  the  fact  that 


n - 1 « 

1 A ( s ) D , , ( s )e<  z A (s)e-*-0  as  N->®. 

v=N  + l v n-v+1  ,n  + l -V  = N + 1 v - - 


Corollary  5.10:  Denoting  the  maximal  eigenvalue  of  a non-negative 
matrix  C by  sp(C),  we  have  that 


sp[lim  Dn>n  + 1(s)]  < sp[G(s )]  (5.11) 


Proof : By  (5.3),  (5.8)  and  the  strict  positivity  of  G(s),  note 

that  R(s)>>0.  Let  n(s)  be  the  spectral  radius  of 

R(s)=TTm  D _ + i(s),  and  let  ;x(s)>>0  be  a right-eigenvector  of 
n-*-®  * 

R(s)  associated  with  n(s).  Then 


n(s)x(s)  = R(s)x(s) 


< I A (s)Rv(s)x(s)  by  Lemma  5.9 

“ v= 0 v 

= z A (s)nv(s)x(s) 

n v 


implying  that 

A* ( n ( s ) ,s)x(s)  > n(s)x(s)  (5.12) 


where 

oo 

A*(z,s)  = I A ( s ) z v , 0<z<l 

v=0  v 


Let  t(z,s)  be  the  Perron-Frobenlus  eigenvalue  of  A*(z,s).  Now 
(5.12)  implies  that 


r 


! 
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C ( n ( s ) , s ) > n(s)  (5.13) 

Now,  for  fixed  s>0,  a theorem  due  to  J.  F.  C.  Kingman  [3] 
implies  that  the  function  1 o g C ( e ^,s)  Is  convex  and  decreasing 
for  t>0.  Also  G(s)  Is  the  unique  solution  of  equation  (5.<0  whence 
it  follows  that  sp[G(s)]  Is  the  unique  solution  zQ  of  the  equation 

£(z  ,s)*z  , 0<z<l  , 

“ 

so  that  zQ=e  , where 

-tn 

log[c(e  u,s)]=-t0,  tQ>0. 
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- t*  -t* 

Setting  s*e  in  equation  (5.13),  we  obtain  logC(e  ,s)>~t* 

which  upon  consideration  of  the  graph  in  Figure  1,  Implies  that 

t*>tQ.  This  clearly  is  equivalent  to  the  Inequality 

sp[R(s  ) ]<sp[G(s  ) ] , which  we  set  out  to  prove. 

Proof  of  Theorem  5.1:  By  Formulas  (5.8)  and  (5.11),  we  have 

spP  im  Dn  n + 1 (s)]>sp[G(s)]>sp[7Tm  D +1(s)] 
n-*®  ’ n-»®  ’ 

But  by  (5.3), 

sp[l_im  Dn  n + 1 ( s ) ]<sp[TTm  0 ,(5)]. 

n-*-®  n->* 

Thus  we  have 

spPim  Dn  n + 1 (s)]=sp[G(s)]=sp[TTm  Dn  ^ ( s ) ] . 
n * n * 

Since  1 i m D , ( s ) < T i m D „.,(s),  and  both  are  irreducible  non- 

n-*’®  n*-« 

negative  matrices,  their  spectral  radii  can  be  equal  only  if  the 

inequality  is  actually  an  (entrywise)  equality  [1].  This  proves 

the  existence  of  lim  0 , ( s ) . Also  (5.8)  an^  the  fact  that 

n-*»  n,n  1 

sp[G(s)]»sp[lJjn  Dn  n+1(s)]  Implies  that  11m  Dn  n+1(s)»G(s). 
n*°°  * * 
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